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16. Abstract A theoretical analysis of the properties of hydrogen has been made to determine 
the suitability of hydrogen as a cryogenic wind tunnel test gas. By using cryogenic 
hydrogen-, instead of air or cryogenic nitrogen, as the wind tunnel test gas, a signifi- 
cant increase in the test Reynolds number may be achieved without increasing the aero- 
dynamic loads. Under sonic conditions, for example, compared to air at ambien-t tempera- 
ture, cryogenic hydrogen at a pressure of one atmosphere produces an increase in 
Reynolds number of a factor of approximately fourteen while cryogenic nitrogen, at the 
same pressure, produces an increase of only a factor of about six. 

The theoretical saturation boundary for parahydrogen is well defined. Thusj any pos- 
sible effects caused by the 1 iquef action. of the test gas can easily be avoided by 
knowing the maximum local Mach number on the model. 

The nondimensional ratios used to describe various flow situations in hydrogen were 
determined and compared with the corresponding ideal diatomic gas ratios. _ The results 
wereused to examine different inviscid flow configurations. This investigation con- 
cluded that the relatively high value of the characteristic rotational temperature 
causes the behavior of hydrogen, under cryogenic conditions, to deviate substantially 
from the behavior of an ideal diatomic gas in the compressible flow regime. Therefore, 
if an ideal diatomic gas, is to be modeled, cryogenic hydrogen is unacceptable as a 
wind tunnel test gas in a compressible flow situation. 

However, at low Mach numbers where the assumption of incompressibility is valid, the 
deviations in the isentropic flow parameters for cryogenic parahydrogen from the 
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Continuation of Item 16. Abstract 

corresponding flow parameters for an ideal diatomic gas are negligible. Thus, in 
the incompressible flow regime, cryogenic hydrogen is an acceptable test gas. 

Hydrogen properties and fan drive-power requirements related to a hydrogen wind 
tunnel were also examined. The drive power requirements were found to decrease 
with decreasing temperature and may be adequately predicted by using modified 
versions of the ideal gas equations. Since gaseous hydrogen is capable of 
penetrating and degrading the mechanical characteristics of numerous materials, 
materials known to be compatible with hydrogen must be used exclusively, in the 
design of a hydrogen wind tunnel to avoid problems as a result of exposure to 
gaseous hydrogen. 

A literature survey resulted in the conclusion that although hydrogen is a highly 
combustible substance, safety codes exist which, when followed, minimize the risk 
involved in handling hydrogen. 
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SmiMARY 

A theoretical analysis of the properties of hydrogen has been 
made to determine the suitability of hydrogen as a cryogenic wind 
tunnel test gas. By using cryogenic hydrogen, instead of air or 
cryogenic nitrogen, as the wind tunnel test gas, a significant 
increase in the test Reynolds number may be achieved without increasing 
the aerodynamic loads. Under sonic conditions, for example, compared 
to air at ambient temperature, cryogenic hydrogen at a pressiure of 
one atmosphere produces an increase in Reynolds number of a factor of 
approximately fourteen while cryogenic nitrogen, at the same pressure, 
produces an increase of only a factor of about six. 

The theoretical sattiration boundary for parahydrogen is well 
defined. Thus, any possible effects caused by the liquefaction of the 
test gas can easily be avoided by knowing the maximum local Mach 
number on the model. 

The nondimensional ratios used to describe various flow situations 
in hydrogen were determined and compared with the corresponding ideal 
diatomic gas ratios. The results were used to examine different 
inviscid flow configurations. This investigation concluded that the 
relatively high value of the characteristic rotational temperature- 
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causes the behavior of hydrogen, under cryogenic conditions, to 
deviate substantially from the behavior of an ideal diatomic gas 
in the compressible flow regime. Therefore, if an ideal diatomic 
gas, is to be modeled, cryogenic hydrogen is unacceptable as a wind 
ttjinnel test gas in a compressible flow situation. 

However, at low Mach numbers where the assumption of incompressibility 
is valid, the deviation in the isentropie flow parameters for 
cryogenic parahydrogen from the corresponding flow parameters for an 
ideal diatomic gas are negligible. Thus, in the incompressible flow 
regime, cryogenic hydrogen is an acceptable test gas. 

Hydrogen properties and fan drive-power requirements related to 
a hydrogen wind tunnel were also examined. The drive power require- 
ments were fovmd to decrease with decreasing temperature and may be 
adequately predicted by using modified versions of the ideal gas 
equations. Since gaseous hydrogen is capable of penetrating and 
degrading the mechanical characteristics of numerous materials , 
materials known to be compatible with hydrogen must be used exclusively 
in the design of a hydrogen wind tunnel to avoid problems as a result 
of exposure to gaseous hydrogen. 

A literature survey resulted in the conclusion that although 
hydrogen is a highly combustible substance, safety codes exist which, 
when followed, minimize the risk involved in handling hydrogen. 
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II'ffiEODUCTIOiv 


The wind tunnel has "been the primary tool for experimental aero- 
dynamic research and development for many decades. Because of the 
increase in size and speed of proposed aircraft configurations, a need 
has developed over the years for improvet] ground testing capabilities 
in terms of Reynolds nxmiber. The need has been well docmented, for 
example, in references 1 and 2, Programs aimed at the development of 
efficient transport aircraft and maneuvering fighter aircraft to operate 
at transonic speeds have demonstrated soma of the deficiencies in the 
testing cajjahilities of present day wind tunnels, especially with 
respect to adequate Reynolds number simulation. In the trarasonic 
region, one of the major prohleras is the inability to adequately 
determine the effect of Reynolds munber on the boundary 3 ayer - shock 
wave interactions, and in turn, on the performance, stability, and 
trim characteristics of the aircraft. 

The test Reynolds number at a given Mach number may be increased 
by using a heavy gas rather than air as the test gas , by increasing 
the size of the model and tunnel, by increasing the operating pressure 
of the tunnel, and by reducing the test temperature. The dinamic 
pressure, mass flow rate, and power consumption of the tunnel will, in 
general, he affected by the method chosen to increase the Reynolds 
number. 

Tne use of a heavy gas, such as Freon-12, is a well-known method 
of achieving high Reynolds numbers. However, when compressibi3.ity 
effects become significant, the difference between rhe ratios of 
specific heats between the heavy gas ana air leads to an improper flow 
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simulation, thus making the use of a heavy gas unsuitable for tests at 
transonic speeds. 

The more common approaches of increasing the size and increasing 
the stagnation pressure in order to increase the Reynolds number involve 
serious problems concerning the cost of construction and operation, the 
model and support loads, and the difficulty of providing continuous- 
flow capability due to the large power reqviirements . 

The fourth method of increasing the test Reynolds number, reducing 
the test temperature, offers an attractive solution to many of the 
preceding problems. However, as the test tonperatiire is reduced the 
properties of the test gas may begin to deviate significantly from 
the properties of air, which, for all practical purposes under the 
conditions encountered in flight, behaves like an ideal diatomic gas. 

Ihe real gas effects may become appreciable at cryogenic temperatiires 
and may alter the aerodynamic test results. Therefore, an analysis 
must be performed to determine if a particular gas at cryogenic 
temperatures and anticipated operating pressures is acceptable for 
wind timnel research. 

A cryogenic wind tunnel concept has evolved at the HASA Langley 
Research Center in which liquid nitrogen is sprayed directly into the 
txinnel circuit to cool the tunnel structure, remove the heat input 
from the drive fan, and balance the heat conducted into the stream 
through the tunnel walls. Nitrogen is the resulting test gas using 
this cooling procedure. Throughout the development of the concept, 
experimental and theoretical studies have been performed to assess the 
suitability of cryogenic nitrogen as a test gas. By analyzing real 
gas isentropic and normal shock solutions , and by measuring the surface 
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presstires on am airfoil, Adcock, Kilgore, and Ray, in reference 3, 
concluded that the real gas effects of nitrogen at cryogenic tempera- 
tures pose no problem in its application in a cryogenic transonic wind 
tunnel operating at stagnation, pressures up to five atmospheres. 

Nitrogen is, of course, not the only gas that may be considered for 
use. at cryogenic temperatures. In theory, any diatomic gas will do. 

One of the most promising diatomic gases in addition to nitrogen is 
hydrogen. Because of the lover vapor temperature, even higher test 
Reynolds numbers may be obtained by using cryogenic hydrogen rather 
than nitrogen. For example, at a pressure of one atmosphere, when 
compared to ambient air, cryogenic nitrogen and hydrogen, produce 
increases in Reynolds number by factors of approximately six and four- 
teen respectively for sonic flows. Alternatively, in a given size 
wind tunnel, it is possible to achieve a desired test Reynolds number 
at a lower total pressure by using hydrogen instead of nitrogen, 
thereby reducing the model, sting, and balance loads. 

This thesis contains a theoretical analysis of those properties 
of hydrogen needed to determine if, in fact, it is suitable for use as 
a wind tunnel test gas. The thesis is divided as follows: In Chapter 

I, some of the fundamental properties essential to understand hydrogen 
are presented. In Chapters II, III, and IV some of the nondimensional 
ratios used to describe various flow situations in hydrogen are 
determined and compared with the corresponding ideal diatomic gas 
ratios. Chapters V and VI present inviscid flow examples to illustrate 
the deviations studied in Chapters II through IV. Hydrogen properties 
and fan drive— power requirements related to a hydrogen wind tunnel are 
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'presented in Chapter VII. Finally, reconmendations are made concerning 
the use of hydrogen as a possible wind tunnel test gas. 
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CHAPTER I 

PROPERTIES OF GASEOUS HYDROGEN 

1.1, Introduction . 

Hydrogen is a homo nuclear diatomic molecule having a relatively 
high characteristic rotational temperature. Because of its molecular 
composition hydrogen demonstrates various anomalies in its properties 
affecting its suitability as a wind tunnel test gas. Tlie purpose 
of this chapter is to introduce the reader to the major peculiarities. 

1.2. Ortho-, Para- a Normal- ^ and Equilibrium Hydrogen . 

Hydrogen consists of a mixtiore of two different types of molecules 
having different optical and thermal properties. These two distinct 
forms of hydrogen are known as parahydrogen (P-I 12 ) and orthohydrogen 
(o-H^)* The parahydrogen molecules have antiparallel nuclear spins 
and even rotational q,uantum numbers, while the orthohydrogen molecules 
possess parallel nuclear spins and odd rotational q^uantum numbers. 
Figure 1 illustrates the difference between ortho- and parahydrogen. 

The thermodynamical eq^ui librium between ortho- and parahydrogen 
is governed by Boltzmann’s distribution law as presented in reference 
4. According to this law the fraction of the molecules, , of the 
total number of molecules, N^, in the rotational state j is given by 
the equation 



where 
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Pj = statistical weight of the state j 
Ej = statistical energy of the state j 
kg = Boltzmann's constant 
T = temperature 


and where may he taken as equal to one. Since all parahydrogen 

molecules have even rotational quantum numbers and all orthohydrogen 
moleciiles have odd rotational quantiom manhers, the ratio, B, of the 
concentration of the two forms is given hy 


e = 




( 1 . 2 ) 


Or, considering the energy of the state j and the statistical weight 
as governed hy the quantum statistics applicable to the nuclei 


E. = i(j+l) = j(j+l)h 
^ dni 


= W.Sj 


23+1 

3(2j+l) 


for j 
for j 


even 

odd 


where 

I = moment of inertia 
h = Planck's constant 

g = statistical weight of a degenerate state 


then 



T 



vhere 

9^ = - 8^1.837 K (from reference h) 

since 

h = 6.55 X 10 erg/sec 
I = U .67 X 10 "^^ g-cm^ 
kg = 1.373 X 10“^^ erg/K 

Eq^uation l.U gives the composition of equilibrivu!» hydrogen as a 
function of temperature and is illustrated in figure 2. At high 
temperatures (around room temperature) the constant h is much less 
than kgT causing the equilibrium ratio to approach the ratio 


3 1/3 (high temperature limit) 


(1.5) 
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I’he high temperature equilibrivun mixture, consisting of twenty- five 
percent pafahydrogen molecules and seventy-five percent orthohydrogen 
molecules , is classified as normal hydrogen and is the equilibrium 
composition of hydrogen at room temperature and above. 

She value of b becomes much greater than as the temperature 

approaches absolute zero, and from equation 1.4 the equilibrium ratio 
approaches the limit 


3 CO (low temperatiire limit) (1.6) 

Uius, at very low temperatures 'equilibrium hydrogen consists of 
essentially all parahydrogen molecules. Hydrogen at the lower 
temperature limit consists of pure parahydrogen because as the 
temperatiire is lowered all molecules gradually pass into the lowest 
state having the quantum number zero, which is a parahydrogen state. 

However, according to reference 4, the rate of conversion from 
parahydrogen to orthohydrogen, that is, the rate at which equilibrium 
between the two species is approached, is extremely slow unless the hydrogen, 
is heated in the presence of a suitable catalyst, such' as oxygen or charcoal. 
In general the self conversion rate (no catalyst, present ) is a function of 
temperature. For example at liquid air temperatxires the half life, 
defined as the time required to convert one half of the excess ortho- 
or para composition present at the starting time, of the conversion is 
greater than a year while at 923 K and 0.0657 atmospheres the half life 
is only on the order of ten minutes. Since the hydrogen that would 
be used as a wind tunnel test gas would probably be stored in liquid 
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form, and considering *that the experimental run times are relatively 
short compared to the para— ortho conversion time, the simplifying 
assumption vill be made that the hydrogen gas being used in the wind 
tunnel is pure parahydrogen. 

1.3 Characteristic Rotational Temperature ^ 

As previously mentioned, the characteristic rotational temperature 
for hydrogen is relatively high. Because of this high value, hydrogen 
can exist in the gaseous state for temperatures lower than the 
characteristic rotational temperature. A range of influence associated 
with the characteristic rotational temperature can be expected so that 
at the lower temperature limit hydrogen should behave like an ideal 
monotomic gas and at the higher temperature limit (around room 
temperature) hydrogen should behave like an ideal diatomic gas. 

The range of influence may be analyzed by considering the rotational 
partition function for a homonuclear diatomic molecule given as (see, 
for example, reference 5) 


Q. 




'rot (odd) 


S (2j- 


2 (2j+l)exp (-o(j+l) 

j=even ( ^ j 


E 


Sn(s +1) 

p (2j+l)exp 

j=odd 


or 

Q. 




rot (even) 


(2j+l)exp <-o(j+l) 


E 

j=even 


(2j+l)exp 


J=odd 


|-j(j+l) 

|-j(j+i) 


(l.Ta) 


(1.7b) 
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where is the nuclear degeneracy. Equations 1.7a and l.Tt are for 

odd and even mass number homonuclear molecules, respectively. 

In most cases the characteristic rotational temperature, 0^, is 
near absolute zero (for example, 0^ is 2.1 K for oxygen and 2.9 K for 
nitrogen). Thus, the temperatures at which the gaseous phase exists 
are much greater than the characteristic rotational temperature. Under 
these conditions equations 1.7a and 1.7b both reduce to the form 



i" Z^(2j+l)exp 
0=0 



( 1 . 8 ) 


The sum for 6^ « T is essentially equal to the area under a curve 
of (2j+l)exp{-j(o+l)0^/T} versus j and the sximmation may be replaced 
with an integral 



CO 

Jo 


(2j+l)exp 



(1.9) 


from which 



= 1 2 T 

2 ^n e 


( 1 . 10 ) 


Because the value of the characteristic rotational temperature for 
hydrogen is large (9^ = 84.837 K for hydrogen), equation 1.10 does not 
hold throughout the gaseous region of hydrogen. Equation 1.7a must be 
used to determined the rotational partition function of hydrogen because 
it has an odd mass number (equal to one). The nuclear degeneracy is 


given by 
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= 2m^ + 1 


{ 1 . 11 ) 


where is the nuclear spin quantum number ^ for hydrogen). 

Substituting into equation l.Ta yields 


Q = Q +0 

rot rot, par a rot, ortho 


j=even 


(2j+l)exp <-j(j+l) 


E 

j=oad 


(2j-M)exp — > (1*12) 


where the first summation represents the contribution from the para- 
states and the second summation represents the contribution from, the 
ortho-states. 

Figure 3 illustrates the influence range of the characteristic 

rotational temperature for parahydrogen. At temperatures well below 

the characteristic rotational temperature the rotational partition 

function, Q , is approximately equal to one, signifying that the 

ro X ,para 

rotational mode is unexcited and that the gas should behave similar 
to an ideal monatomic gas. For parahydrogen this occurs for tempera- 
tures below approximately 75 K. For temperatures well above the 
characteristic rotational temperature the rotational partition function 
should approach the ideal diatomic gas rotational partition function 
given by equation 1.10. For parahydrogen, the rotationeJ. partition 
function must be multiplied by four since at higher temperatures 
equilibrium hydrogen consists of only twenty-five percent parahydrogen. 
Thus, for temperatures above approximately 200 K, the rott.tional 
contribution is the dominant factor and parahydrogen may be treated 
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as an ideal diatomic gas. Therefore, parahydrogen may he expected to 
behave similar to an ideal monatomic gas for temperatures less than 
75 K and similar to an ideal diatomic gas for temperatures greater than 
200 K. 

The main part of this thesis will be concerned with how the 
characteristic rotational temperature influences the properties of 
gaseous hydrogen and the resulting deviations of hydrogen from the 
behavior of an ideal diatomic gas. 

1.4 Parahydrogen Equation of State. 

There are literally hundreds of "equations of state" wliich have 
been developed to describe parahydrogen. For example, references 6, 7> 
8, 9s and 10 describe some of the equations of state to be found in the 
literature. 

Today there are four general methods developed by the National 
Bureau of Standards (UBS) used to computerize the thermodynamic 
properties of fluids. These are; linear interpolation; thirty-three 
term modified Benedict Webb Rubin (MBWR); polynomial interpolation; 
and a twenty term I'ffiWR. The linear interpolation method is an 
interpolation on a table of property values stored in computer. The 
method is fast but not very accurate (average of 1 %) and requires more 
computer core storage than the other methods. The most accurate 
method, better than 0.1^ average, is the polynomial interpolation 
method. The major disadvantage of the polynomial method is' that the 
computations are slow relative to the other methods. 

A good compromise, and a widely used method, is the thirty-three 
term modified Benedict Webb Rubin equation of state. The method is 
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faster than the polynomial interpolation and has an accm^acy only 
slightly vorse { 0 . 1 % average )• The major disadvantage of the is 

the great inaccuracies near the critical region ( 20 % error in density 
and 5^ in temperature). 

On a recommendation from the personnel of the NBS Cryogenics 
Division, the polynomial interpolation equation of state for para- 
hydrogen vas selected for use in this study and a computer program vas 
obtained directly from them. The prograia gives the thermodynamic and 
related properties of parahydrogen from the triple point to 300 K at 
pressures to 1000 bar. 

1.5 Saturation Boundary . 

The minimum total temperature vhich could conceivably be used in a 
hydrogen vind tunnel depends on several factors, but is basically 
dependent upon the onset of condensation. The assumption is made that 
condensation must be avoided under all conditions to be encountered 
during the testing of a model in the vind tunnel and that condensation 
is most likely to occur in localized lov pressure regions near the 
model. These lov pressure regions may be thought of in terms of high 
local Mach numbers encountered near the model. To calciolate the 
theoretical saturation boundary, the assumption is' made that the test 
gas expands i sent Topically from the total conditions to the maximum 
local Mach number. 

Using the vapor pressure equation for parahydrogen given in the 
NBS program, along with the above assumptions, the theoretical 
minimum total temperatxires have been calculated for a range of totad 
pressures and maximum local Mach numbers, and are presented in 
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figure 4. For example, if a hydrogen wind tunnel was operating at a 
total pressure of five atmospheres and a maximum local Mach number of 
0.5* the minimum total temperature which could he used to avoid 
condensation would be approximately 28 K. However, if the maximum 
local Mach number was increased to 1.5, the total temperature would 
have to be increased to approximately 37 K. 

Total conditions may be selected to avoid any possible errors in 
the data caused by condensation by using figure 4 and by knowing the 
maximum local Mach number expected dxxring a test. The theoretical 
total conditions necessary to avoid any possible errors due to 
condensation are believed to be conservative since the test gas remains 
at the high local Mach numbers for only very short periods of time. 

If the localized high Mach number regions are small, a nonequilibrium 
frozen flow would exist where the test gas is supersaturated. 

As previously mentioned, the Eeynolds niimber will increase with 
decreasing temperature. Thus the maximimi test Reynolds number which 
can be achieved at a given total pressure is determined by the 
saturation boundary. This will be discussed in greater detail at the 
end of this chapter. 

1.6 Thermal and Caloric Imperfections . 

The need to consider the effects of thermal and caloric imper- 
fections on flow characteristics for parahydrogen at cryogenic 
temperatures is apparent from figures 5 and 6 where the compressibility 
factor, Z, and the ratio of specific heats, C^/C^, as obtained from 
reference 11, are shown over a range of temperatures. Departure of Z 
from the ideal gas value of unity is designated as a thermal imperfection 
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vhereas a departure of C /C frora tlie ideal diatoihic gas value of 

Jr ^ 

1.^ is desigjiated as a caloric imperfection. 

Figure 6 illustrates hov the specific heat ratio is affected by the 
characteristic rotational temperature. For temperatures above 200 K 
the specific heat ratio for parahydrogen is approximately the sair.e as 
the ideal diatomic gas value of l.lf whereas j, for tempex'aturos less 
than 75 K the specific heat ratio corresponds more closely to the ideal 
monatomic gas value of 1,67. 

To evaluate the significance of the nonideal beliaA^ior> real-gas 
isentropic and real-gas normal-shock flow solutions, along with other 
flow situations , have been obtained for parahydrogen for a range of 
pressures and temperatures down to saturation. The calculation 
procedure and discussion of the results for these flow situations are 
presented in the following chapters. 

I . 7 Transport Properties . 

The two transport properties that are of most concern in 
experimental aerodynamic research are the viscosity, n, and the thcnnal 
conductivity, k. These properties are relatively independent of small 
changes in pressure and are plotted in figures 7 ^i^d 8 as a function 
of temperature as obtained from reference 11, 

The two properties along witn the specific heat at constant 
pressure combine to give the Prandtl number, Pr, an important dimen- 
sionless parameter in viscous heat transfer problems, defined as 
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This parameter for parahydrogen is shown as a function of temperature 
at different pressures in figure 9 as obtained from reference 11. 

Eucken's relation, derivable through kinetic theory (see, for 
example, reference 12), gives the Prandtl number as a function of 
the specific heat ratio, S'S 


Pr 


W 

9y - 5 


(l.lU) 


The relationship gives a value of 0.667 for the Prandtl number of an 
ideal monatomic gas where y = 1.67 and a value of 0.737, which is 
approximately the correct value of 0.725, for the Prandtl number of ' 
an ideal diatomic gas where y = l.U. 

Comparing Eucken's relation to figure 9 and considering the 
characteristic rotational temperature influence, there appears to be a 
major inconsistency. One explanation for the deviation is that 
Eucken's relation is a reasonably good approximation only for gases at 
ordinary temperatures. At cryogenic temperatures, the transport 
properties do not vary in the same manner as they do at room temperatures 
thereby causing Eucken's relation not to hold at cryogenic temperatures. 

Cryogenic temperatures influence the transport properties in such 
a way as to cause the Prandtl number at one atmosphere and around 
25 K to be equal to the ideal diatomic gas value of 0.725 instead of the 
expected value of 0.667, the ideal monatomic gas value. So, under 
these conditions, the thermal boundary layer profile may be expected 
to vary from the velocity boundary layer profile for parahydrogen in 
the same manner as an ideal diatomic gas. 



IT 


As noted in the introduction, it is necessary to match, as nearly 
as possible, the Reynolds numher between flight conditions and the 
conditions maintained during aerodynamic testing in. the wind tiuinel. 
With the ever increasing flight Reynolds mmibers, much work has been 
devoted to increasing the test Reynolds number by decreasing the test- 
ing temperature, by increasing the testing pressure, and by using 
different gases (see, for example, reference 13). Figure 10 
illustrates the benefits, from the viewpoint of increased Reynolds 
nmber per meter, of using cryogenic hydrogen compared to nitrogen or 
air assuming an isentropic expansion from the total conditions. For 
example, at a Mach number of 0.1, the Reynolds number per meter for 
hydrogen at the total conditions of one atmosphere and 25 is about 
fifteen times that of air at standard conditions and approximately 
three times that of nitrogen at the total conditions of one atmosphere 
and 100 K. 

The Reynolds number per meter for hydrogen and nitrogen are 
compared at a free stream Mach number of 1.0 and a total pressure of 
one atmosphere at various total temperatures in figure 11. The figure 
illustrates the large benefits that may be obtained by using hydrogen 
instead of nitrogen from a Reynolds mmiher standpoint. The benefits 
are a direct result of the lower total temperatures which can he used 
before the onset of condensation occurs. 

Thus, the maximum obtainable Reynolds number per meter is directly 
dependent upon the saturation limits as imposed by the maximum local 
Mach number generated on the model being tested. This is illustrated 
for hydrogen in figure 12 where the Reynolds number per meter is shown 
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as a function of total temperature for various total pressures with 
the theoretical saturation boundaiy superimposed as a function of the 
maximum local Mach number, The figiure shows that if the maximxim 

local Mach number is l.U a factor of 17.5 increase in Beynolds number 
may be obtained at a pressure of five atmospheres by operating at 
cryogenic temperatures rather than ambient temperatures. This repre- 
sents a Reynolds number which is fojrty-five times that of air at 


standard conditions. 
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CHAPTER II 

ISENTROPIC FLOW OF PARAIIYLEOGEN 


II »1 Introduction. 

A useful approximation to portions of flows in aerodynamic re- 
search is the isentropic flow of an ideal diatomic gas. Consequently, 
any candidate test gas should have isentropic flow ratios similar to 
an ideal diatomic gas. The purpose of this chapter is to present the 
real gas isentropic flow solution and results for parahydrogen. 

II. 2 Isentropic Flow Solution . 

By definition an isentropic process occurs under the condition of 
constant entropy, and the solution may he calculated from the stagna- 
tion conditions. An analytical solution may be formulated from the 
simple equation of state for an ideal gas. The equation of state for 
a real gas is more complex than the state equation of an ideal gas and 
an analytical solution for the real gas is not as easily determined. 

As noted in Chapter I the equation of state for parahydrogen used 
in the isentropic flow solution was developed by the Rational Bureau 
of Standards. The computer program for the equation of state was 
modified for the flow solution into a subprogram form (identified as 
subprogram THERMO) where the inputs are the temperature and pressure, 
and the outputs are the thermodynamic and related properties. 

In addition to the main program, a subprogram (identified as 
ISEHT) was also developed that makes use of THERMO. ISEWT calculates 
the pressure and other thermodynamic properties for a given 
temperature and entropy by using a modified interval halving technique, 
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and it is this subprogram that the main isentropic flov solution 
program centers around. 

The main program iterates on the Mach nxaaber and employs an 
interval halving technique, modified for the isentropic flow solution, 
to calculate the isentropic expansion of parabydrogen at various 
total temperatures and total pressures. Generalized flov charts 
describing the main program and subprogram ISEira are presented in 
figure 13 and the steps involved in the solution are discussed in 
detailed in reference l^t. 

II. 3 Isentropic Expansion Coefficient . 

The pressure and density relationship for an isentropic expansion 
of an ideal gas is described by the exponential equation 

p “ (2.1) 

where y is the ratio of specific heats and is a constant for an 
ideal gas. As noted earlier in figure 6, the specific heat ratio 
for the real gas parabydrogen is not constant but with decreasing 
temperature deviates considerably from the ideal diatomic gas value of 
l.li. Because the specific heat ratio varies with temperature, 
equation 2.1, with y equal to a constant, would not be expected to 
be valid for the isentropic flow solution of parabydrogen. 

In reference 15, Woolley and Benedict indicate that equation 2.1 
may hold true for an isentropic expansion of a real gas except that 
the exponent would no longer be equal to the specific heat ratio. They 
define the exponent as the isentropic expansion coefficient, 5, which 
may be calculated from the formula 
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Z 


1 


T7 



( 2 . 2 ) 


and indicate that equation 2.1 may he used for an isentropic expansion 
of a real gas provided that ^ is used instead of y* 

In this thesis the isentropic expansion coefficient was determined 
from the expression 



where the states i-1 and i represent an increment along an isen- 
trope which is equivalent to a variation of 0.05 in Mach number. The 
variation of the isentropic expansion coefficient with total 
temperature at a total pressure of one atmosphere is illustrated in ' 
figure lU for expansions of parahydrogen to a Mach number of 2.0. This 
figure illustrates that the isentropic expansion coefficient is 
approximately constant outside the characteristic rotational temperature 
influence range and equivalent to the specific heat ratio for an ideal 
monatomic and an ideal diatomic gas below and above the influence range 
respectively. These values agree with those predicted by Woolley and 
Benedict when the rotational degrees of freedom are either unexcited 
or fvilly excited. 

II. Isentropic -Flow Ratios . 

As noted in section II. 2, since, by definition, an ideal diatomic 
gas is both thermally and calorically perfect , the isentropic flow 
ratios may be easily solved for by using the simple ideal gas equation 
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of state. These fundamental equations are found to he 



1 + 0.5(y - 1)M^ 
1 + 0.5(y - 
1 + 0.5(y - 1)M^ I 





(2.3) 

(2.U) 

(2.5) 

( 2 . 6 ) 


and are a function of Mach number, M, and the specific heat ratio, Y> 
only. 

As previously noted, any candidate wind tunnel test gas should have 
behavior similar to air which, under flight conditions, acts for all 
practical purposes like an ideal diatomic gas. For this reason, the 
isentropic flow ratios for parahydrogen will be compared with the ideal 
diatomic gas ratios deterained from equations 2.3 to 2.6 with the spe- 
cific heat ratio, y» equal to 1.4. 

The deviation of the isentropic flow ratios from the gas values 
are presented in graphical form. The plots were obtained from the 
tabulated values generated by the computer program discussed in 
reference l4. 

The deviation of the pressure ratio for an isentropic expansion 
of parahydrogen to Mach 1.0 from the ideal diatomic gas value is 
illustrated in figure 15. The figure shows that the deviation is 
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mainly dependent upon the total temperature variation while the effect 
of the total pressiire on the deviation is comparatively small. Once 
again the real gas behavior of parahydrogen, this time in terms of the 
isentropic expansion pressure ratio, may he explained by the influence 
of the relatively high value of the characteristic rotational tempera- 
ture since the greatest deviations occur within and below the 
characteristic rotational temperature influence range, defined in 
section 1.3 as 75 K to 200 K. 

The effect of Mach number on the deviation of the isentropic 
pressure ratio is shown in figure l6. The deviation is shown to be 
insignificant for low values of Mach number, which is expected because 
the flow may be assvimed to be incompressible luider these conditions. 

It is interesting to note that the maximum deviation of the isentropic 
pressure ratio for parahydrogen occurs at approximately the sairie Mach 
number as that determined for nitrogen by Adcock in reference l6. This 
occurance surely must be coincidental because Adcock shows that for 
nitrogen this deviation is in the positive sense and the maximum value 
increases with both Mach number and pressure whereas for parahydrogen 
the deviation is in the negative sense and the maximum value increases 
with pressiare but decreases with Mach number. Figure l6 also verifies 
the fact that the deviation with respect to the total pressure varia- 
tion is compeuratively small. 

Figure 17 shows how the temperature ratio for an isentropic 
expansion of parahydrogen to Mach 1.0 deviates from the ideal diatonic 
gas value. Like the pressure ratio, the deviation in temperature 
ratio is seen to be only slightly dependent upon the variation in total 
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pressure and to be basically dependent upon the total temperature. 

The isentropic temperature ratio as a function of total temperature 
curve is identical in shape to the isentropic pressmre ratio versus 
total temperature curve. The influence of the characteristic rota- 
tional temperature is again the dominant factor responsible for the 
deviation. 

The Mach number effect on the isentropic temperature ratio 
deviation is illustrated in figure l8. Similar to the isentropic 
pressure ratio, the deviation of the temperature ratio is insignificant 
in the incompressible range of lov Macb number values. This figure 
also verifies that the deviation with respect to the total pressure is 
negligible when compared with the deviation caused by the total 
temperature variation. 

An illustration of the isentropic density ratio deviation for an 
expansion of parahydrogen to Mach 1.0 is given in figure Ip. This 
figure shows that the density ratio deviates in the same manner as 
the isentropic pressure and temperature ratios. That is, the isentropic 
density ratio varies mainly with respect to the total temperature and 
only comparatively slightly with respect to total pressure. The 
characteristic rotational temperature influence is once again the 
dominant factor controlling the deviation. 

The deviation of the isentropic density ratio is shown to be 
generally less than the deviation of the pressure and temperature 
ratio. The deviation may be explained by realizing that pressure, 
density, and temperature are all interelated. As previously noted, 
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the isentropic pressiire ratio deviation and the isentropic temperature 
ratio deviation have curves vhich are similar in shape, so that vhen the 
density ratio deviation is calcxzlated, the deviations of temperature 
ratio and pressure ratio have a tendency to cancel each other, thereby 
keeping the density ratio deviation closer to unity. 

Figure 20 gives the effect of Mach number on the deviation of the 
isentropic density ratio. In the incompressible range of low Mach 
numbers the real gas effects are once again shown to be insignificant- 
Figure 20 also verifies that the deviation with respect to the total 
pressure variation is comparatively insignificant when considering 
the deviation caused by the total temperature variation. 

The real gas effects on the isentropic stream-tube area ratio for 
parahydrogen are illustrated in figures 21 and 22. Since A -is 
equivalent to when the Mach number is equal to 1.0, the deviation 

from the ideal diatomic gas value is zero at sonic conditions and 
small in the transonic region, as sho\m in figiire 21. The deviation 
is shown to be more dependent on total temperature than on total 
pressures and is smaller in the incompressible range than in the 
supersonic range. The deviation with respect to total temperature, 
illustrated in figure 22, is again dominated by the characteristic 
rotational temperature influence. 

II , 5 Conclusions . 

1. The isentropic expansion coefficient for parahydrogen cannot 
be assumed to be constant but is governed by the 
characteristic rotational temperature influence range. 
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2, The characteristic rotational temperature is also found to he 
the dominant factor controlling the deviations of the 
isentropic flow parameters* 

3. The deviations of the isentropic flow parameters are primarily 
a function of the total temperature, with the influence of 
the total pressure being small in comparison. 

4* At low Hach numbers, where the assumption of incompressibility 
is valid, the isentropic flow deviations are negligible* 
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CHAPTER III 

PRAI'TDTL-^MEYER EXPAHSIOH FLOW OF PARAHYDROG-EH 
XII * 1 Introduction > 

One type of isentropic flow which requires further discussion is 
the supersonic isentropic flow that occurs in the turning of a flow 
around a corner* This type of flow is known as a Prandtl— Meyer flow* 
The relation between the flow inclination and change in Mach number in 
an isentropic turn is given by the Prandtl-Meyer function. 3y knowing 
the flow inclination all of the isentropic flow ratios may be obtained 
at the resulting Mach number through the use of the Prandtl-Meyer 
function* The purpose of this chapter is to present the Prandtl— Meyer 
flow solution and the dependency of Mach number on the Prandtl-Meyer 
function for flow in parahydrogen as compared to the results for an 
ideal diatomic gas* 

III *2 Prandtl-Meyer Expansion Flow Solution, 

Consider the semi-infinite flow field above a convex wall as shown 

in figure 23 where the subscripts a and b are used to identlf^'* the 

upstream and downstream conditions of the expansion wave. Here a 

uniform parallel supersonic flow of Mach M is expanded to a higher 

a 

speed, by an increase in the local flow angle, 9^ - 9^. As in all 
Prandtl-Meyer flows, the Mach lines, which are Inclined relative to 
the local flow direction at the Mach angle, p, given by 

- 1/2 


tan p = (M^ - 1) 


(3.1) 
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form as straight lines from the wall. In addition, the flow parameters 
are constant along these lines. Therefore, the problem of determining 
the flow above the waJ.1 transforms into the determination of the 
velocity and the thermodynamic state as a function of the initial 
conditions and the change in the local flow direction imposed by the 
wall. 

The governing differential equation for Prandtl-Meyer flow can be 
determined from the geometry of a differential velocity change across 
a Hach line, considered as a weak, discrete expansion wave, as shown 
in figure 2h, Here the tangential components of the velocities on 
both sides of the Mach line are equal, that is, 

u^ = (u + 5u)^ - (3.2) 


From the law of sines 


u + 6u 
u 



or, using trigonometry 


6u _ cos(u) 

u cos{ii)cos(dO) - sin(vi)sin(60) 


(3.3) 


(3.U) 


So, for small 60 


1 + ^ = cos(y) 

u cos(y) - sin(p)d0 


(3.5) 


This equation reduces to 
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1 + ^ = 


u 1 - tan(ji)d6 


(3.6) 


or, neglecting higher order terms 


1 + “ = 1 + tan(p)d6 


(3.7) 


Combining equation 3*7 with equation 3.1, the well-known differential 
equation for Prandtl-Meyer flow is obtained 


d8 = - 1 ^ 

u 


(3.8) 


To obtain a solution, two additional equations are required, the 
adiabatic energy equation, 


1 2 


1 2 




(3.9) 


and the thermodynamic equation of state, given here in terms of the 
speed of sound, c, for an isentropic process, 


c = c(H,S) = c(K,S„) (3.10) 

a 

With the initial conditions given, equations 3.8, 3.9, and 3.10 can 
be solved for any three of the four variables, 0, u, c, H, in terms of 
the fourth. 

For an ideal gas, equation 3.10 can be written in a simple form 
and the equations can easily be solved. For a real gas, the equation 
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of state is more complex than the ideal gas state equation and such a 
simple solution to the system of equations is not as easily determined. 

A computer program was written to solve the set of equations for 
parahydrogen. The main part of the program follows the generalized 
flow chart given in figure 25 and the steps involved are discussed in 
reference IT* 

III. 3 Prandtl-Meyer Expansion Flow Results . 

As previously stated, the system of equations for an ideal gas 
may easily be solved for, resulting in the well-laiown equation for 
the Prandtl-Meyer function for an ideal diatomic gas 

= 2.kk93 tan"^ (^0. ko825 ~ ^ (3.11) 

The Frandtl-Meyer expansion flow results >for parahydrogen will he 
compared with the ideal diatomic gas values since, as mentioned 
earlier, any candidate wind tunnel test gas must hehave like air which, 
unaer flight conditions, acts like an ideal diatomic gas* 

The deviations from the ideal diatomic gas values for the Prandtl- 
Meyer expansion flow in parahydrogen are presented in graphical form. 

The results, shown in figure 26, were obtained from the tabulated 
values, given in reference 17s generated by the computer program 
mentioned in the previous section. 

Figure 26 shows the Mach number deviation as a function of total 
temperature for Prandtl-Meyer expansions from tlach 1.0 tlirough angles 
of five and twenty-five degrees at total pressures of one and ten atmospheres. 
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As was the case for the isentropic flow ratios, the deviation is mostly 
dependent upon the total temperature variation while the deviation with 
respect to the variation in total pressure is comparatively small. 

Once again the influence of the characteristic rotational temperature 
appears to he the dominant factor affecting the Mach number deviation, 
just as it is in the deviations of the isentropic ratios. The Mach 
n\mher deviation also appears to increase with increasing Prandtl- 
Meyer function, 0. This is explained hy the facts that the Mach number 
increases with 0 and, from Chapter II, that the isentropic deviations 
increase with increasing Mach number. 

III.U Conclusions. 

1. The characteristic rotational temperature is the dominant 
factor influencing the Mach number deviation for the 
Prandtl-Meyer expansion in parahydrogen . 

2. The Mach n\anber deviation is small with respect to total 
pressure variations when compared to the influence of total 
temperature. 

3. The Mach number deviation generally increases with increasing 
values of the Prandtl-Meyer function. 
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CHAPTER IV 

NORMAL SHOCK WAVES IN PAEAHYDROGEN 
IV »1 Introduction. 

Shock waves are usually present around models tested at transonic 
and supersonic speeds and the investigator woiild like the jump condi- 
tions across the shock wave to have a behavior similar to air. So, 
any candidate test gas should have normal shock ratios similar to an 
ideal diatomic gas. The purpose of this chapter is to present the 
normal shock flow solution and results for parahydrogen . 

IV. 2 Normal Shock Wave Solution. 

- Consider a standing normal shock wave as shown in figure 27- If 
the subscripts a and b are used to identify the upstream and 
downstream conditions of the shock wave, the conservation equations of 
mass, moment imi, and energy are 




(U.l) 


P + p M^c^ = P, + p, M?c, 
a a a b '^b b 


(4.2) 


0^9. j D 


(4.3) 


Equation 4.3 may be written in terms of the static enthalpy as 


1 ,.2 2 _ „ ^ 1 ^2 2 

+ T;* M C = H, + 7T M, C, 

a 2aa b 2ob 


(i^.4) 
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Using equation k.l, equations k.2 and may Ue written as 


~ P - P /p-^^) 

D a "^a a a ' a U 


( 1 ^. 5 ) 




a 2 a a 


1 - 


(h.6) 


If the conditions upstream of the shock are known, equations 4.5 and 
4.6 constitute two equations for the three unknowns P^, p^. The 

system of equations can he completed hy assuming the gas to he in 
thermodynamic equilihrium on each side of the shock wave, which will 
supply the equation of state 


P =’p(P,H) (4.T) 

The requirement given hy the second law of thermodynamics , that the 
entropy downstream of the shock wave, S^, he greater than the entropy 
upstream of the shock wave, S^, assures the uniqueness of the solution. 
For an ideal gas, equation 4.7 may he written as 

P = P[t/(y - 1)]/H (4.8) 

where y is the specific heat ratio. The system of equations may then 
easily he solved for analytically. 

In the case of a real gas the equation of state is more complex 
than the state equation of an ideal gas, and an analytical solution of 
the system of equations cannot he determined algebraically. 
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A computer program was written to solve the normal shock wave 
system of equations for paraliydrogen. The equations were solved hy 
iterating on the upstream Mach number and employing a modified 
interval-hadving method to converge on the downstream static enthalpy 
and density. To determine the conditions on either side of the shock 
wave, the program used the assumption of local thermodynamic 
equilibrium and employed subprograms TPiEEMO and ISEKT, which were 
discussed in Chapter II. ' 

A generalized flow chart for the main part of the computer program 
is given in figure 28 and the steps involved are discussed in reference 
l8. One step that should be mentioned here is step 10 where the static 
temperature downstream of the shock is initially assumed. For this 
assumption the ideal, diatomic gas 'equation was utilized, that is 


T . 
a,i 


T 


t ,a 


(1 + 0 . 21 ' 4 ^ 






0.2M^ 

a. 


.-1 




T 




ih.9) 


where the expression enclosed "between the absolute value signs 
represents the deviation from the ideal diatomic gas value for the 
previous iteration. The above method provided rapid convergence 
throughout the domain studied. 


IV. 3 Normal Shock Wave Results. 

The normal shock wave ratios for an ideal gas may, as previously 
noted, be solved for analytically and are given as 


2Y^^ - (y - 1) 

a 


Y + 1 


ih.9) 
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T 

a 




- (Y - 1) 


(y - I)i'r + 2 


(y + 1)^M^ 


^a • (y - 1)M^ + 2 


M 


(Y - l)M^ + 2 

\'[2Yl4j - (Y - 1) 


1/2 


(U.IO) 




(U.12) 


X 


^t,b _ Pt‘,b 

(y + l)M^ 

a 

Y-1 

1 


(Y - l)M^ + 2 

— Si ^ 


1 

> 

1 

CM ^ 


1 

Y-1 


(4.13) 


T 


:^=i 

t,a 


(4.14) 


The fundamental equations are found to be functions of the upstream 
Mach number, M , and the specific heat ratio, y» only. 

As mentioned earlier, any candidate wind tunnel test gas should 
behave similar to an ideal diatomic gas. For this reason, the normal 
shock wave ratios for parahydrogen will be compared with the ideal 
diatomic gas ratios determined from equation 4.9 to 4.l4 with the 
specific heat ratio, Y» equal to 1.4. 

The deviations of the normal shock wave ratios for parahydrogen 
from the ideal diatomic gas values are presented in graphical form. 
The plots were obtained from the tabulated values given in reference 
l8 which were generated by the computer program discussed in 


section IV. 2. 
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The deviations of the static ratios P, /P_, T, /T , and p, /p , 

D Da Da 

across a normal shock wave in parahydrogen from the ideal diatomic gas 
values at various upstream total pressures and for an upstream Kach 
numher of 2*0 are illustrated in figures 29 to 31* The figures show 
that the nonideal gas behavior of parahydrogen is mainly dependent 
upon the upstream total temperature variation, while the deviations with 
respect to the upstream total pressure are comparatively small* The 
greatest influence on the deviations is the previously mentioned high 
value of the characteristic rotational temperat-ure of hydrogen* 

Figures 32 through 3^ show how the deviations of the static ratios 
across a normal shock wave in parahydrogen vary \rith respect to the 
upstream* Mach number, M^, at various upstream total temperatures and 
for an upstream total pressure of one atmosphere. The figures show 
that the deviations generally increase with increasing upstream Mach 
number. They also show that the deviations, as a function of total 
t empe rat Tire , have the same shape at the various upstream Mach nimibers, 
and appear to be greatly influenced by the relatively high value of the 
characteristic rotational temperature* ^ 

The real gas effects on the Mach number downstream of a noriaal 
shock in parahydrogen are illustrated in figure 35 for an upstream 
Mach number of 2.0* The deviation from the ideal diatomic .gas value 
is shown to be only slightly dependent upon the upstream total pressure 
when compared to the dependence upon the upstream total temperature. 

Once again the dependency is explained by the characteristic rotational 
temperature range of influence. 
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Figures 36 through 38 give the 'effect of the upstream total 

pressure on the deviations of the total ratios s , /T, , 

tjh t^a t>D t,a 

P+ -u/P+ > from the ideal diatomic gas values for an upstream Mach 

u ) u b ^ a 

number of 2.0. The figures show that the total pressure ratio and the 
total density ratio deviations are mainly dependent upon the upstream 
total temperature variation and are only comparatively slightly 
influenced by the variation in the upstream total pressure. The 
deviations are shown to be predominantly influenced by the 
characteristic rotational temperature. 

On the other hand, the deviation of the total temperatvire ratio 
from that of an ideal diatomic gas does not appear to be influenced 
by the characteristic rotational temperature. The greatest deviations 
of the total temperature is on the same order of magnitude as the total 
pressure and density deviations. 

Figures 39 through 4l illustrate how the deviations of the total 
ratios across a normal shock wave in parahydrogen vary with respect 
to upstream Mach n\miber,,M , at various upstream total temperatures 

8i 

and for an upstream total pressure of one atmosphere. The figures 
show that the deviations generally increase with increasing upstream 
Mach number. Comparing figure 40 with figure 3Ts the deviation from 
the ideal diatomic gas value of the total temperature ratio is 
shown to be only slightly dependent upon the upstream Mach nxuaber 
when compared to the dependency of the upstream total temperature. 

IV. 4 Conclusions. 

1. The characteristic rotational temperature is the dominant 
factor controlling the deviations of the static ratios 
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acrpss a normal shock wave in parahydrogen from the ideal 
diatomic gas values. 

2. The deviations of the normal shock wave static ratios increase 
with increasing upstream Mach number* 

3. The deviation's of the static ratios are small with respect to 
the total pressure variations when compared to the influence 
of the total temperature. 

U. The characteristic rotational temperature is also the dominant 
factor controlling the deviations of the total pressure and 
total density ratios across a normal shock wave in parahydrogen 
from the values for an ideal diatomic gas^ while the total 
temperature ratio deviation is not influenced by the 
characteristic rotational temperature. 

5* The deviations of the total pressure and total density ratios 
are small with respect to the total pressure variations when 
compared to the influence of the tojbal temperature, while the 
inverse is true for the total temperature ratio. 

6. The deviations of the stagnation ratios across normal shock 
waves increase with increasing upstre^. Mach number. 
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CHAPTER V 

NOZZLE FLOW OF PARAHYDROGEN 


Y*1 Introduction ♦ 

The previous chapters have described hov hydrogen differs from an 
ideal diatomic gas for various flov situations. Also, » the dramatic 
increase in Reynolds number attainable through, the use of cryogenic 
hydrogen instead of air or cryogenic nitrogen was shown in Chapter II, 
along with the conclusion that cryogenic hydrogen will behave like an 
ideal diatomic gas in the incompressible flow regime. 

This chapter considers inviscid one-dimensional cryogenic 
parahydrogen flow, both isentropic and with a normal shock, through 
the diverging section of a supersonic nozzle. Through this analysis 
the differences between cryogenic hydrogen and an ideal diatomic gas 
will be studied in specific compressible flow situations. 

V.2 Isentropic Flow. 

The dimensions for the diverging section of a Mach 2.0 nozzle 
were obtained from reference 19. The values were then fitted with an 
equation of the form 



where 

L total length of the diverging section of the nozzle 
X distance downstream from the nozzle throat 
A area at distance x 


area at the nozzle throat 



•uo 


. and the values of are as follows: 

7 q = l-O 

y^ = 0. 0169041 

Yg = 6.873320 
Y3 = --18.309T 
v\ = 23.9338 
= -16.5135 

yg = 4.6773if 

The area distribution ^ along with the resulting one-dimensional isen- 
tropic flow ratios for an ideal diatomic gas, is presented in figure h 2 . 

The computer program discussed in reference l4 was modified to 
obtain the one-dimensional, inviscid, isentropic flow of parahydrogen 
through the diverging section of the -Mach 2.0 nozzle given by equation 
5-1. Since the purpose of this study is to investigate the behavior of 
hydrogen xmder cryogenic conditions where the Reynolds number advantage 
occurs, the isentropic flow ratios were calculated for total conditions 
of 1^5 K and one ahKiosphere . These conditions give a Reynolds nuuiber 
per laeter value of ahout. seven times that of air at one atmosphere and 
300 K. 

The results at the above total conditions are given in figure 43 
in tenas of the deviations from the ideal diatomic gas values. As 
expected from the resialts of Chapter II, the deviations increase with 
increasing x/L and show a rather large deviation in the isentropic 
ratios even at the throat. Tlie figure leads to the conclusion that 
although substantial increases in Reynolds numbers may be obtained with 
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cryogenic hydrogen, hydrogen is nnacceptahle as a test gas if it is 
essential to model an ideal diatomic gas iinder supersonic isentropic 
flow situations. 

V.3 Mozzle Flow vith a IIoiTnal Shock 

The isentropic nozzle flov solution discussed in the preceding 
section vas altered to study the nonideal gas behavior of parahydrogen 
associated with a one- dimensional normal shock wave. This was done by 
using the real gas normal shock wave relations derived in Chapter IV 
and assuming isentropic flow upstream and downstream of the shock. As 
in the previous section, to talte advantage of hydrogen *s high Reynolds 
number capability at cryogenic temperatures, the Mach number distribu- 
tion and nozzle flow ratios were calculated for upstream total 
conditions of 45 K and one atmosphere. 

The parahydrogen and ideal diatomic gas results for shocks . 

occTirring at Mach 1.3 and 1.8 are shown in figures 44 and 45, 
respectively. As expected from the results of Chapter II, the deviation 
in the shock wave location between parahydrogen and an ideal diatomic 
gas increases with increasing shock Mach number. This is predictable 
from the fact that in the supersonic range, the deviation for para- 
hydrogen from the ideal diatomic gas value of the isentropic stream- 
tube area ratio, A*/A, increases with increasing Mach number (as shown 
in figure 22). The results indicate that the shock wave location on an 
airfoil being tested at supersonic speeds in cryogenic paraliydrogen 
would probably be in error by an unacceptable amount, when compared 
to the shock wave location in ,air, thereby possibly nullifying any 
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advantage of testing at the increased Reynolds numher which may be 
obtained by using cryogenic par ahydrogen . 

The results of figures 44 and 4^ are in agreement with the results 
of Chapter II, and demonstrate that the deviations in the Mach number 
and nozzle flow ratios are generally smaller downstream of the shock 
wave than they are upstream of the shock wave* The results support 
the conclusion reached in the Chapter II that the isentrcpic flow 
ratios for cryogenic parahydrogen at low subsonic Mach numbers do not 
deviate substantially from the ideal diatomic gas values. 

V.4 Conclusions. 

1. Cryogenic hydrogen is unacceptable as a test gas if it is 
essential to model an ideal diatomic gas under supersonic 
isentropic flow situations. 

2. Cryogenic hydrogen is unacceptable as a test gas if it is 
essential to model an ideal diatomic gas under supersonic 
flow situations where a shock occurs. 



CHAPTER VI 


FLOW ABOUT A DIA140KD-SHAPED AIRFOIL 
VI ♦! Introduction . 

Chapter V concluded that, when compared to an ideal diatomic gas, 
cryogenic parahydrogen is unacceptable as a supersonic wind tunnel test 
gas. The purpose of this chapter is to verify this conclusion and, 
like the nozzle flow study, to gain insight into the flow characteris- 
tics of cryogenic parahydrogen in the transonic Mach n\im‘ber region. 

This is done hy using the shock-expansion theory to study the inviscid 
flow about a diamond-shaped airfoil. 

VI. 2 The Oblique Shock Wave . 

One of the basic concepts involved in the shock-expansion theory 
is the oblique shock wave, the geometry of which is given in figure 46. 
The oblique shock wave may be treated in the same manner as a normal 
shock wave taking into consideration the additional velocity component, 
V, tangent to the shock wave. 

The continuity equation across an oblique shock may be given in 
terms of the nomal velocity components as 



where u is the velocity component normal to the shock wave and p is 
the density. 

From the oblique shock wave geometry given in figure 46 the 


following may be written 



tan (3 - 6) = 


( 6 . 2 ) 


tan 3 = 


(6.3) 


Using tlie trigonometric identity 


tan(3 - 6) = .-rf 

1 + tan3 tan6 


equations 6.2 and 6.3 may be combined to yield 


1 tan6 

^2 u 

— = ( 6 , 

u, u ' 

1 1 

1 + — tan6 

V 

But, the tangential component of the velocity, v, may be written as 


/ 2 2 
V = - u^ 


so that 


Introducing the speed of sound, c, the above may be written 


= /M? 


"l J ^Vl 


1 


. ( 6 . 8 ) 



where is the upstream Mach n'umber and the ratio referred 

to as the equivalent normal shock Mach number* Substituting equation 
6.8 into equation 6.5 gives the following result based only on the 
geometric considerations of an oblique shock wave 



From the above equation and the equation of continuity, 6.1, the 
solution for the oblique shock wave may be determined at a specified 
value of 6. 

VI. 3 The Shock-Expansion Theory . 

The oblique shock wave along with the isentropic expansion wave 
are the two fundamental components used to analyze many two-dimensional 
flow problems by uniting the appropriate combination of the two 
solutions. The shock- expan si on theory may be used to analyze the flow 
over two-dimensional airfoils. 

Consider, for example, the symmetrical diamond- sect ion airfoil in 
figiire Uj . A shock forms at the leading edge compressing the flow to 
pressures and above and below the airfoil. Centered 

expansion waves located at the shoulders expand the flow to pressures 
^3 ^5* FliiQ-llys the trailing edge shock recompresses the flow 

to (nearly) the free-stream pressure value. 

A supersonic wave drag is induced on the airfoil due to the over- 
pressure on the forward faces and the underpressure on the rearward 
faces. This type of supersonic drag exists even in the idealized. 
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nonviscous fluid and is fundament ally different from the frictional 
drag and separation drag associated -with the viscous boundary layer. 

For unit span^ the supersonic wave drag associated with the diamond- 
section airfoil is given hy 

= [(Pg + P]^) - (P3 + P5)]t (6.10) 

where t is the maximm thickness of the airfoil. 

VI . 4 Flow Solution ♦ 

Now consider the symmetric diamond-section airfoil in figure 4 T 
with a chord length, £, of three meters flying at an altitude of 
8000 meters and a Mach number , of 1 * 3 . At this altitude, reference 
20 gives the following conditions 

static pressure: 0.352 atm 

kinematic viscosity: 2.904 x 10 m /sec 

speed of sound: 308.07 m/sec 

static density: O.526 kg/m^ 

From the above values, the Reynolds umber based on the chord length 
is 4 l .4 X 10 ^. To match the Reynolds number in a cryogenic hydrogen 
wind tunnel, with a total pressure of one atmosphere and a total 
temperature of 50 K a chord length of 0.4576 meters is required, 
which could easily be accommodated in a wind tunnel. 

The lift and drag coefficients associated with the diamond-section 
airfoil given in figure 47 may be determined using the shock-expansion 
theory. For the ideal diatomic gas the values are easily determined 



using the equations given in chapters II, III, and IV following the 
calciilatlon procedure outlined in reference 21. 

Because of the complex equation of state for a real gas , an 
analytical solution for hydrogen is -not as easily determined. The 
results obtained in chapters II, III, and IV may be used along with the 
shock-expansion theory to calculate the lift and drag coefficients. 

Using the normal shock wave solution tables, given in reference l8, 
along with equations 6.1 and 6.9»-the Mach number normal to the shock 
wave upstream of the shock wave may be determined for various deflection 
angles, 6. Figure 48 gives the solution for the free-stream conditions 
of one atmosphere and 50 K, and for deflection angles, 6, of two and 
four degrees, and an upstream Mach number, M^, of 1.3. 

The lift and drag coefficients were calculated for various free- 
stream Mach numbers , , airfoil half angles , (}) , and angle of attack , 

a, and are presented in table I. 

VI. 5 Comparison of Results . 

The errors in the lift and drag coefficients do not vary in any 
consistent manner. Occasionally the errors are less than one percent, 
but, in general, the errors are large enough to make hydrogen 
unacceptable as a cryogenic wind tiinnel test gas if an ideal diatomic 
gas is to be simulated. 

VI . 6 Conclusions . 

The results from the diamond-shaped airfoil study confirm the 
conclusions reached in the nozzle flow study. That is, considering 



that high local Mach number regions may occur in transonic and super- 
sonic flow situations, and recognizing the possibility of large errors 
in the determination of the lift and drag coefficients, cryogenic 
hydrogen is unacceptable as a transonic and supersonic ' wind tunnel 
test gas if an ideal gas is to- be simulated. 
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CHAPTER VII 

SELECTED WIND TUNl^EL DESIGN CONSIDERATIONS 
VII >1 Introduction , 

The purpose of this chapter is to discuss selected design aspects 
vhich must he considered if a cryogenic vind tunnel is to he built or if 
an existing wind tunnel is to he modified to use cryogenic parahydrogen 
as the test gas* No attempt is. made to cover all of the design 
considerations associated with a cryogenic wind tunnel. Only those 
aspects peculiar to the use of hydrogen will he examined. These 
include: drive power reguirements , cooling requirements, material 

compat ability^ and safety. 

VII. 2 Drive Power Requirements . 

This section consists of a theoretical analysis of the power 
required for isentropic compressions of cryogenic parahydrogen to 
determine how the real gas effects influence the drive power require- 
ments of a closed circuit wind tunnel. The analysis covers a 
temperature range of 25 K to 300 a pressure range of one to five 
atmospheres, and a fan pressure ratio range of 1.019 to 1.200. The 
values cover the range of variables, likely to he encountered in a fan- 
driven cryogenic wind tunnel with hydrogen as the test gas. The results 
are compared to air at total conditions of one atmosphere and 300 K. 

In reference 22, Adcock and Oghum present power calculations 
for the isentropic compressions of cryogenic nitrogen, and their 
analytical wind tunnel model, presented in figure 49, shall he 
examined for the hydrogen case. 
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The assumption is made that the only losses that occur are between 
the tunnel throat and fan. The simplifying assumption is a direct 
result of wind tunnel data which show that most of the tunnel losses 
do occur in this portion of the tunnel because of the higher flow 
velocities. As a consequence of the simplification, the total condi- 
tions are the same at the fan outlet and the t\mnel throat . 

To remove the heat conducted through the walls of the tunnel and 
the heat added to the stream by the fan, the assumption is made that 
a cooling system is placed upstream of the fan. A further assumption 
is made that the cooling is accomplished by injecting liquid hydrogen 
directly into the t\mnel circuit thereby cooling the circuit through 
evaporation. The location of the cooler is upstream of the fan rather 
than downstream to allow a more thorough mixing of the evaporating 
hydrogen with the main stream. In the case of nitrogen, Adcock and 
Ogbum have shown that the additional mass flow due to the liquid nitrogen 
being injected for cooling is at most two percent of the mass flow. 

Their studies indicated that while the absolute power levels were 
increased by about two percent because of the additional mass being 
compressed, the ratio of the power required to compress nitrogen to 
that required to compress an ideal diatomic gas was insignificantly 
affected. Since they were primarily interested in this ratio, the 
assumption was made that cooling occured without mass addition. The 
same assumption will also be made here, realizing that the power 
calculations do not include the additional power required to compress 
the mass flow of the coolant. 

Based on the fan pressure ratios necessary to achieve a given test 
section Mach number in several existing tunnels, the fan pressure 
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ratios presented in table II have been assumed for this analytical 
study. 

Based on experience TTith existing tunnels ^ the fan pressure ratio 
needed to achieve a given test section Mach nimiber is assumed to be 
invariant vith stagnation temperattire and pressure. In accordance 
vith the conclusions made so far vith respect to the usefullness of 
hydrogen as a test gas, the low subsonic Mach number region will be 
studied in greater detail than the higher Mach number region. 

The tunnel mass flow calcixlations are carried out at the throat 
conditions. For subsonic velocities the assumption is made that the 
throat and test section Mach numbers are identical. For supersonic 
speeds, the effective flow area in the test section must be greater 
•than that of the throat. In practice, this larger effective area is 
Created either by diverging the walls of the test section or by 
allowing some of the test gas entering the test section to flow through 
slotted or porous test section walls into a plenum chamber surrounding 
the test section. In the latter case, the mass may be removed from 
the plenum chamber by auxiliary suction or it may be allowed to reenter 
the tunnel circuit at the diffuser entrance. For this analysis 
it will be assuifted that all the mass that passes through the tunnel 
throat will also pass through the fan. 

The test gas of a closed circuit fan driven wind tunnel is forced 
to flow around the txmnel circuit by the energy imparted by the fan. 
Assuming that the steady flow compression that occurs at the fan is an 
isentropic process, the power per unit area imparted to the test gas is 
given by 
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where 




m = mass flow rate per unit area of the test gas 


(H^ 2 - = enthalpy difference, at constant entropy 

* ’ s 

across the fan. 

The mass flow rate per unit area may he determined hy 


m = pMc 


and for an isentropic process of an ideal gas, this may he expressed 
as a function of the stagnation conditions and Mach number 




Y-^-l 

2{1-y) 


Also, for an ideal gas the enthalpy difference may he expressed as 


(H - II ) = T. ^ I 

t,2 t,l ^ p t,l[\P^^^ 


where (Pt,2/^t,l^ is coEimonly referred to as the fan pressure ratio. 
Equations T»3 and 7^^ may then be substituted into 7*1 to give the 
pover requirements for an ideal gas. 

In the case of a real gas such a closed form solution cannot be 
as easily determined. But, by using the isentropic flov computer 
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program discussed in chapter II ^ and specifying the fan pressure ratio, 
the throat Hach nurnher , and the total conditions at the throat , the 
power req^uirements for a wind tunnel using parahydrogen as the test 
gas may be calculated using the previously discussed analytical model. 

Figure 50 gives the power requirements for a hydrogen wind tunnel 
relative to the power requirements for an ambient air wind tunnel as 
a function of fan outlet total temperature, T . Each fan outlet total 
pressure line, corresponds to every test section Mach number at 

which the power requirements were calculated. The maximum deviation 
was approximately one-half of one percent between a Mach niunber of 1.2 
and a Mach number of 0.1. This result suggests that the power require- 
ment ratio is relatively independent of Mach number, a fact which 
shall be applied later. 

Table III shows how the test section size, the tunnel pressure, 
and the tunnel power requirement is affected by changing the test gas. 
The table illustrates the advantage of hydrogen compared to nitrogen 
or air* At a given Mach number and Reynolds number the drive power 
requirement for hydrogen is less than the power requirement for nitrogen 
or air. The advantage of hydrogen is a direct result of the saturation 
boundary allowing operation at a greatly reduced total temperature. 

Since the advantage of parahydrogen lies in its high Reynolds 
number capability, the power requirements for a certain Reynolds number 
advantage will be considered. The power requirement is analyzed best 
as a ratio of the parahydrogen power requirement to the power required 
for air near standard conditions. Kilgore, in reference 13, suggests 
that both the Reynolds number and power are directly proportional to 
the total pressure. Therefore, if the Reynolds number, divided by the 



total pressure, is plotted against the required power, divided by the 
total pressure, only one curve should arise since the total pressure 
dependency should cancel out • Figure 51 gives this function for 
parahydrogen in terms of the Reynolds number advantage (ratio of the 
Reynolds number for parahydrogen at various total pressures and 
temperatures to the Reynolds number for air at one atmosphere and 
300 K), divided by the total pressure, as a function of the required 
power ratio (power per unit area for parahydrogen at various total 
pressures and temperatures divided by the power per unit area for 
air at one atmosphere and 300 K), divided by the total pressure. The 
values for air were, determined by using the ideal gas equations for 
an ideal diatomic gas. 

The solid line labeled ideal gas in figure 51 evolved from a 
modification of the closed form solution of the ideal gas equations. 
The Reynolds number per meter is given by 


(Re /meter) = 


pMc 

n 


(7.5) 


By using the isentropic flow equations for an ideal gas this equation 
may be given by 


(Re /meter) = 


n 



1 + ^ 
2 


„ -X_ 
Y-l 


(7.6) 


The viscosity, n» is usually given by Sutherland's equation, but, as 
pointed out in reference 7, a more accurate value for parahydrogen is 
given oy Biller's foi*mula (see the appendix). The static temperature, 
T,' is given by 
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(7.7) 


As noted in chapter I, the specific heat ratio, y, for parahydrogen 
is not constant but varies with temperature. For this analysis the 
specific heat ratio is given as 



k=0 


(7.8) 


where 

T = temperature, K 

and where the values of Y, are as follows; 

k 


^0 = 

2.6735 


^1 = 

-1.8562 

X 

H 

0 

1 

ro 

^2 = 

6.9001 X 

0 

H 

^3 = 

-1.0288 

X 

10“ 5 

^4 = 

5.3581 X 

10 ^ 

Y5 = 

-6.6509 

X 

10-12 

II 

VD 

-8.744 X 

10”13 

^7 = 

2.894 X 

10-15 

^8 = 

-2.9464 X 10"^® 






Equation 7.8 approximates the specific heat ratio, C^/C^, given in 
figure 6. 

The drive power requirements for parahydrogen were then approxi- 
mated by using equations 7.1, 7.3, and 7.4 with y = 1.4 and 

the appropriate values for the specific heat at constant pressure, C , 

F 

and the gas constant, R. The cixrve was calculated for a Mach number 
equal to one. Calciilations were performed at other Mach numbers but 
these resulted in similar curves demonstrating the relative independence 


of Mach number. 
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VII > 3 Cooling Requirements ♦ 

The amoimt of liquid hydrogen required for cooling may "be divided 

* 

into tvo categories: i) the amoimt initially required to cool the 

cryogenic timnel to its operating temperature^ and, ii) the amount 
required to remove the heat being conducted through the tunnel vails 
and the heat energy added to the stream by the drive fan vhile the 
tunnel is running* 

The amount of liquid hydrogen required to cool a cryogenic wind 
tunnel to its operating temperature is dependent upon the cool-down 
procedure used as well as the physical characteristics of the tunnel. 

To calculate the minimum liqxiid requirement, the simplifying 
assumptions are made that no heat is conducted through the insulated 
tunnel walls and that no heat is added by the drive fan to the test gas 
during the cool-down process. With these assumptions the minimum 
liquid hydrogen requirement for cool-down occurs when the tunnel is 
cooled slowly in such a way that the hydrogen gas will leave the tunnel 
circuit at the same temperature as that of the warmest part of the 
tunnel structure. Assuming further that all of the liquid hydrogen 
vaporizes, the maximum amount of liquid hydrogen required for cool-down 
may be calculated when only the refrigeration available in the latent 
heat of vaporization is utilized and the hydrogen gas leaves the tunnel 
at the saturation temperature. 

If a represents the mass of liquid hydrogen required to cool a 

unit mass of material through a given temperature range, then, 

following the analysis given in reference 23% O . = 0.031 and 

min 

^ ” 0.l4 for cooling a stainless steel structure from 260 K to 


30 K. 
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Cooling through such a large temperature range would not always 
occur since \inder some circumstances the tunnel circuit would be 
allowed to remain cold between runs. Under these conditions the heat 
to be removed would be equal to the heat gained by the tunnel 
structure eind test gas by conduction through the insulation. This heat 
conduction could be kept to an acceptable level by properly designing 
the insulation. 

One possible way of reducing the cost of cooling a cryogenic 
hydrogen wind tunnel would be to use relatively inexpensive liquid 
nitrogen to cool the structiore to approximately 77 K and to use hydrogen 
to cool the structure to the lower operating temperatxire- 

The heat to be removed while the tunnel is running consists of the 
heat conduction through the walls of the tunnel and the heat energy 
added by the drive fan. The cooling thermal capacity of parahydrogen 
is equal to the latent heat of vaporization (445-5 joules per gram at 
the normal boiling point) together with the sensible heat of the gas 
phase between the saturation temperature and the desired total 
temperature. 

Figure 52 gives the cooling thenaal capacity for parahydrogen for 
a range of final gas conditions. Here the cooling thermal capacity 
is equal to the latent heat together with the sensible heat of the gas 
phase between the saturation temperatvire and the final gas total 
temperature. 

Kilgore, in reference 13 » gives the cooling thermal capacity for 
nitrogen. By comparison, the cooling thermal capacity for 
parahydrogen is much greater than that of nitrogen, on a mass basis. 
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At a temperature of 300 K, the cooling capacity of parahydrogen is 
approximately nine times that of nitrogen and^ on a mass basis, the 
latent heat of vaporization of parahydrogen is more than tvice that 
of nitrogen. 

The additional power required for cooling is best illustrated by 
an exaiaple. The mass flow rate of the liquid hydrogen (LHg) required 
to remove the heat of compression is equal to the drive fan power added 
to the stream, Pwr, given by equation 7,1, divided by the cooling 
capacity, given in figure 52 and Appendix I. 

If e represents the energy required to produce a unit mass of 
liquid hydrogen, then the power equivalent of liquid hydrogen used to 
remove the heat of compression is 

where is the mass flow rate of the liquid hydrogen. — 

Table IV compares the results for cryogenic hydrogen, cryogenic 
nitrogen, and air at standard temperature at a free stream Mach number 
of 1.0, a Reynolds num'ber of 50 x 10^, a total pressure of 2.5 
atmospheres and a local saturation I&ch num'ber of 1.2. The value of £ 
for nitrogen vas obtained from the work reported hy Kilgore in 
reference 13. The value of £ for hydrogen was then selected to be 
consistent with the liquefaction work requirements for nitrogen and 
hydrogen given in reference 2k, The results show that although the 
drive power requirements for cryogenic hydrogen is su'bstantially less 
than the drive power requirements for cryogenic nitrogen or air at 
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ambient temperature, the total power requirements for cryogenic 
hydrogen is considerably greater because of the liquefaction power 
requirements • 

VII, 4 Material Compatibility. 

Atomic hydrogen is capable of entering steel and many other metals 
and alloys# When this penetration occurs, any one of seyeral undesirable 
phenomena may take place which are collectively referred to as hydrogen 
embrittlement* Hydrogen embrittlement is a delayed— failure phenomenon 
which is also known as ** hydrogen-induced, delayed brittle failure" or 
hydrogen-stress cracking" since a loss of ductility as measured in the 
tensile test is not necessarily associated with the condition* 

The presence of hydrogen affects the mechanical behavior of iron 
and steel principally by: reducing the ductility (embrittlement), 

lowering the fracture stress, and causing a delayed brittle failure 
under suitable conditions* If an apparatus is to be used in a hydrogen 
environment it should be constructed of material that would not be 
affected by the hydrogen* 

As a general rule for the, selection of metal alloys to use in a 
hydrogen environment, the aluminum alloys and the stable austenitic 
stainless steels, such as SAE 321, have been found to be insensitive to 
the presence of gaseous hydrogen and, therefore, any system using these 
materials should be free from degradation due to hydrogen embrittlement* 
Thus, if these or similar hydrogen compatible materials are used 
exclusively in a hydrogen wind tunnel, problems caused by exposure to 
gaseous hydrogen will not occur* 
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VII. 5 Hydrogen Safety . 

Perhaps the greatest safety hazard encountered when handling 
hydrogen is its high comhustahility. The ordinary burning of a mixture 
of hydrogen and air can occvir within the broad limits of about four 
percent to seventy-five percent of hydrogen by volume. When combustion 
occurs rapidly an explosion results. Thus unless suitable precautions 
are taicen, a flammable mixture is possible in a hydrogen wind ttinnel 
during some stage of operation and any source of ignition could then 
cause serious trouble. Because hydrogen is a highly corabustable 
substance, the wind tunnel designer and technician should be fajuiliar. 
with the combustion of hydrogen. Tiie following discussion serves only 
as an introduction to this subject and in the event that a hydrogen 
wind tunnel is constructed, the designer should refer to specific 
details in other sources, for example, reference 25. 

V/hen hydrogen and air are burned the reaction is 

^ 2 +^ 0 ^ + I.89W2 = 1*89^2 ^ (7.10) 

This is the stoichiometric mixture of thirty percent hydrogen by volume 
which when burned uses up all of the oxygen in the air. The heat of 
combustion converts the water to steam and raises the temperature of 
the resulting mixture. 'frJhen this occurs in a closed volume, the 
temperature increases from T^^ to T^ and the number of moles 
changes froio n^^ to n^ causing a pressxare change of 



n T 

2 2 

n T 

1 1 


(7.11) 
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Wie above equation assumes an adiabatic process, uniform temperature 
distribution throughout the mixture, and negligible dissociation of 
the combustion products , resulting in a higher calculated pressure 
than actually occurs. 

It is recommended that any wind tunnel designed to use hydrogen 
as the test gas be constructed so that if the hydrogen is to be 
discharged into the atmosphere, it could be burned vmder .controlled 
conditions. The principle is the same as that used in oil refineries 
where combustable gases are burned rather than allowed to accumulate. 

Several safety codes, such as references 26, 27, and 28, exist 
which could be of help during the design of a hydrogen wind tunnel and 
auxiliary equipment. Generally, the building in which the wind tunnel 
is contained should have: l) a vent system for the apparatus, 2) ade- 

quate ventilation, 3) an adequate electrical groxxnding system, and 
no pockets near the ceiling which covild trap hydrogen. 
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CliAPTER VIII 
CONCLUSIONS 

A theoretical analysis to determine the suitability of hydrogen 
as a cryogenic wind ttinnel test gas has been made. The major 
conclusions to be drawn from the analysis are as follows: 

1. A wind tunnel using cryogenic hydrogen j instead of air or 
cryogenic nitrogen, as the test gas will have a significant 
increase in test Reynolds* number without increasing the 
aerodynajnic loads. 

2. The theoretical- saturation boundary for parahydrogen is well 
defined. Therefore, any possible effects caused by the 
liquefaction of the test gas can easily be avoided provided 
that the maximum local I'i'Iach number on the model is known. 

3* The relatively high value of the characteristic rotational 
temperature causes the behavior of hydrogen, under cryogenic 
conditions, to differ substantially from the behavior of an 
ideal diatomic gas when considering the compressible flow 
regime. Therefore, if it is essential to model an ideal 
diatomic gas, cryogenic hydrogen is unacceptable as a wind 
tunnel test gas in compressible flow situations. Consequently, 
a cryogenic hydrogen wind tunnel would not provide an adequate 
simulation of transonic and supersonic flow situations in 


ambient air. 
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4. At low Mach numbers, where the assumption of incompressibility 
is valid, the deviations in the isentropic flow parameters for 
cryogenic parahydrogen from the ideal diatomic gas values are 
negligible. Thus, in the incompressible flow regime, 
cryogenic hydrogen is an acceptable test gas. Consequently, 
a cryogenic hydrogen wind tunnel would provide an adequate 
simulation of incompressible flow situations in ambient air. 

5* A cryogenic hydrogen wind tiinnel requires more power to 
operate than a cryogenic nitrogen wind tunnel. Although a 
smaller size cryogenic hydrogen wind tunnel may be designed 
because of the Reynolds number advantage of hydrogen, more 
power is required to produce liquefied hydrogen than is 
required to produce liquefied nitrogen. This difference in 
the liquefaction power requirements results in the above 
conclusion . 

6. In the event that a wind tunnel is constructed which would use 
hydrogen as a test gas, materials compatible with hydrogen 
must be used exclusively to avoid possible problems as a 
result of exposure to gaseous hydrogen. 

7 . Although hydrogen is a highly combustable gas , safety codes 
exist which, when followed, minimize the risk involved in 
handling hydrogen. 
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Figvire 2. Composition of equilibrium hydrogen. 




Figure 3. Characteristic rotational temperature influence range as illustrated 
Djr the rotational partition function. 
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90 X 10"7 

80 

7r 

60 

VISCOSITY, 

T], 50 

kg/m-sec 

40 

30 

20 

I 

ir 


Figure 






Figure 0- Thermal conductivity, k, for parahydrogen. 
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Figure 11. Reynolds number per meter comparison between hydrogen and nitrogen at 
a free streaia Mach number of 1.0 and a total pressure of one atmosphere. 




Figure 12. Reynolds number restrictions due to theoretical saturation boundaries 
at a free streain Mach number of 1.0. 



































Figure 13. Concluded. 
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Figure 15. Pressure ratios for the isentropic eiicpansion of parahydrogen to a Mach 
number of 1.0 relative to the ideal diatomic gas value. 
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Figure l6. Pressure ratios for the isentropic expansion of parahydrogen to 
various Mach nunhers at a total temperature of 4o K, relative to the ideal 
diatomic gas value. 
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Figure l8* Temperature ratios for the isentropic expansion of parahydrogen to 
various Mach numbers, relative to the ideal diatomic gas value* 




Figure 19 . Density ratios for the isentropic expansion of parahydrogen to a Mach 
nuBiber of 1.0 relative to the ideal diatoBiic gas value. 








Figure 22. Isentropic stream-tube area ratio for parahydrogen relative to the 
ideal diatomic gas value. 




Figure 23. Prandtl-Ileyer flov over a convex vail. 
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Figure 25. Computer program flow chart for the 
Prandtl-Meyer flow solution. 








Figur 8 26 • Mach number deviation as a function of botal ’temperature for Prandtl— 
Meyer expansions from a Mach number of 1.0 “through angles of 5 Q-ud 25 degrees in 
parahydrogen* 
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Figure 27 


Standing normal shock vave. 
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Figure 28, Computer program flow 
normal shock wave solution. 
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Figure 28 . Concluded,. 
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Figure 30. Static temperature ratio across a normal shock wave in parahydrogen 
relative to an ideal diatomic gas for an upstream Mach number of 2.0. 




gure 32. Static pressiire ratio across a normal shock vave in parahydrogen 
Latiye to an ideal diatomic gas for an upstream total pressure of one atmosphere* 
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Figure 35* Mach number downstream of a n.ornal shock wave in parahydrogen relative 
to an ideal diatomic gas for an upstream Mach number of 2.0, 
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x'igure 39* Total pressure ratio across a normal shock wave in parahydrogen relative 
to an ideal diatomic gas for an upstream total pressure of one atmosphere. 
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Figiire k2. Concluded. 
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Figure 43. Deviations of parahydrogen from the ideal diatomic gas for the one- 
dimensional, inviscid, isentropic flow through the diverging section of a Mach 2.0 
nozzle at total conditions of 45 K and one atmosphere. 
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Figure Mach nmber distribution and nozzle flow ratios 

for the one-dimensional, inviscid flow through the diverg- 
ing section of a nozzle with a normal shock occurring at a 
Mach number of 1.3 at upstream total conditions of 45 K 
and one atmosphere. 
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Figiire 45. Mach number distribution and nozzle flov ratios 
for 'the one— dimensional ^ inviscid flow through the diverg- 
ing section of a nozzle with a normal shock occurring at a 
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Figure 46. Oblique shock vave geometry. 



Figure U7. Diamond airfoil geometry. 
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e 48. Oblique shock wave solution for the hydrogen case. 
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Figure 49* AnaJ.y1:ical model of a wind tunnel. 
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Figure 50. Hydrogen wind tunnel power requirements as a function of fan outlet 
total temperature. 



Figure 51. Hydrogen wind tunnel power requirements 
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Figure 52. Heat absorbing capacity of liquid parahydrogen. 
Liquid initial at a pressxire of one atmosphere and a 
temperature of 20.268 K. 
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APPEIIDIX B 
Tables 

TABLE PAGE 

I. Diamond- shaped airfoil results 132 

II. Fan pressure ratio necessary to achieve a given test 

section Mach nvimber 133 

III. Wind tunnel requirements for various test gases at 

free stream saturation. . 13^ 

IV, Wind tunnel requirements for various test gases at 

a free stream Mach number of 1.0, a Reynolds number 

6 / 

of 50 X 10 , a total pressure of 2; 5 atmospheres, 

and a local saturation Mach number of 1.2 135 
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TABLE I. 

Diamond-shaped airfoil results. 
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TABLE II. 

Fan pressure ratio necessary to achieve 


a given test section Mach number 
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TABLE III. 

Wind "tunnel requiremen'ts for "various "tes"t gases a"t free s"fcreani satura"tion 
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TABLE IV. 


Wind tunnel requirements for various test gases at a free stream Mach number 
of 1.0, a Reynolds niunber of 50 x 10^, a total pressure of 2.5 atmospheres. 


and a local saturation Mach nxmber of 1.2. 
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APPENDIX C 

Fixed Point Properties of Parahydro’gen 

Below is a summary, taJcen from reference 11, of the P-p-T data at 
selected fixed points for parahydrogen. 

Critical Point 
T = 32.976 + 0.05 K 
P = 12.759 atm 
P = 31.43 kg/m^ 

Nonnal' Boiling Point 

T = 20.268 K 

P = 1 atm 

p( liquid) = 70.78 kg/m^ 
p( vapor) = 1.338 kg/m^ 

Triple Point 
T = -13.803 K 
P = 0.0695 atm 
p(solid) = 86.50 kg/m^ 
p (liquid) = 77.03 kg/m^ 
p( vapor) = 0.126 kg/m^ 
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APPENDIX D 
Virial Coefficients 

One equation of state for parahydrogen is given in terms of the 
virial coefficients. The virial coefficients are commonly defined as 

P = ETp[l + B(T)p + o(T)p^ + ...] (D.l) 

where B(t) and C(t) are virial coefficients of a power series 

expansion in density. Equation D.l adequately describes the P-p-T 
surface for densities up to about one-half of critical. Reference 11 
gives B(t) and C(T) as: 

For temperatures below 100 K 

B(t) = (b^T + bg + bg/T + bj^/T^)/RT (D.2) 

where 

\ ~ I.93977J+I X 10^ bg = -1.9279522 X 10^ 

b^ = -2.2890051 X 10^ bj^ = I.IO9U088 X 10*^ 

For T < 55 K 

C(T,) = (Cj^T^ + CgT + c^ + Cj^/T + c^/T^ + Cg/T^)/RT (D.3) 

where there are two sets of the c^ 'coefficients. One for temperatures 
below T^ = 32.95 K 
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= 1.051H7T6 X 10^ 

” I.0431U11 X 10 ^ 

= 5. 1405848 X 10^^ 


Cg = -1. 6597141 X lo"^ 


cj^ = -3.2538718 X 10 
cg = -3.3123453 X 10' 


10 


and a second set for temperatures between and 55 K, 

= 1.6971294 X 10^ Cg = -5.0854223 x 10^ 


- 6.7284118 X lo”^ 


cjj = -3.8045171 X 10" 


= 1.0789413 X 10^^ cg = -1.1515642 X 10 


12 


For 55 < T < 100 


where 






( /^ \ 
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- exp 

."■3 

r ' w 

®1 = 

388.682 


^2 

= 45.5 

^3 = 

-0.6 
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= 20.0 

= 

4.0 





a 


5 


(D.4) 


For temperattires above 100 K 


b(t) = 

i=l ^ ^ 


(•2i-4)/4 


(D.5) 


and 


C(T) - 


0 

o 2 


!-5 - 


exp ( 1 - x" 


(d.6) 


where 
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= k2.^Q\ 

= -2.2982 
= IO9.T8I/T 

= 1310.5 


\ = -37.1172 
\ = -3.0484 
Xg = 20.615/T 
= 2.1486 


For equations D.2 through D.6 the units of temperature, T, is Kelvin, 
the units of B(T) are cm^/mole and the tinits of C(T) are 
(cm^/mole)^. 

The uncertainty of B(t) is a maximum of about five percent at 
the highest and lowest temperatures. The uncertainty of C(T) is a 
minimum of around five percent between 55 and 100 K and as much as 

approximately twenty percent for temperatures below the critical 
temperature. 



APPENDIX E 

Latent Heat of Vaporization 


The latent heat of vaporization is the amount of heat required to 
convert a unit mass of a substance from the liquid state to the vapor 
state at a constant pressure and may easily he obtained hy subtracting 
the liquid enthalpy from the enthalpy of the vapor at the same 
temperatTore and pressure. From reference 11, the latent heat of 
vaporization for parahydrogen ranges from Ij 48.2 Joules/gram at the 
triple point to Ult5.5 Joules/gr^ at the normal boiling point to 0 
Joules/grara at the critical point. The uncertainty in these values 
is estimated to be 1,2 Joules/gram. 
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APPENDIX F 

Thermal Conductivity of Parahydrogen 

The thermal conductivity coefficient relates the transfer of heat 
through a material via molecular interaction caused ty a temperature 
gradient across the material. That is 

q = -k grad T (F.l) 

where q is the heat flux, grad T is the temperature gradient and k 
is the thermal conductivity coefficient. 

One of the most popular methods of correlating thermal conductivity, 
for example reference 11, is to separate the property into a mmiher of 
additive parts 


k = k^(T) + kj^(p,T) + k^(p,T) (F.2) 

where is the dilute gas contribution (a function of temperature 

only, kj^(p,T) is the excess function, and k^(p,T) gives the 
enhancement due to the critical point behavior. 

The thermal conductivity of a gas may also be expressed in terms 
of its specific heat at constant volume, C^, its viscosity coefficient, 
n, and its molecular weight, by 


k = k-nC /M 
1 V w 


(F.3) 
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where the factor allows for the fact that the mean energy of the 

molecToles appears larger as the molecules possessing large amoiints of 
kinetic energy travel faster and thus transfer their energy more 
rapidly. T-Then the molecules possess only translational kinetic 
energy (i.e. in the case of hydrogen at temperatures less than 50 K) 
kj^ is equal to 2.5* If the rotationsQ. kinetic energy is completely 
available (i.e. for hydrogen above 300 K) then k^ is equal to 1.9. 
In reference 4, Farkas gives an approximation of k^ as 

2.25 R + C 

c ^ 

V 

where R is the gas constant. Therefore ^ the thermal conductivity 
coefficient is given as 

(2.25 R + c^)n 

M 
w 


k = 


(F.5) 
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APPENDIX G 


Viscosity of Parahydrogen 


The viscosity of a fluid may he separated into two additive parts 


p. = n^(T) + rig(p,T) 


(G,l) 


where Pq(T) is the dilute gas contribution and Pg(p,T) is the 
excess or dense gas contribution and is a strong function of density. 
In the lover temperature range of T < 100 K, reference 7 gives 


r\^{T) = a^^ + a^) (T + a2)/(T + aj^) 


(G.2) 


= 8.5558 X 10' 

= 650.39 


ag = 19.55 
a,^ = 1175.9 


and for the excess viscosity contribution 


n„(p>T) - A(p) exp B(p)/T 


(G.3) 


(Aj_p + Agp + A^p^) 

(1.0 + AjjP + A^p^ + Agp^) 


X 10 


A(p) = 


(G.4) 



1U4 


= 306.U636 


Ag = -3350.628 


and 


B = 


A^ = 38680.92 

Aj^ = -I8.U783 

A^ = 110.915 

5 

Ag = 25.3524 

r \ "bo' 


b ■ 

(p/Bg) ^ - (p/B ) ^ 

+ B^ exp 

Bj^(p/ 3 g) ^ 

1 




(G.5) 


Bq = 10.0 

B2 = 0.07 
= -58.75 

bg = 1.5 


B, = 


= 


b, = 


b^ = 


7.2 

-17.63 

6.0 

3.0 


For equations G.2 through G.5, the temperature, T, is in Kelvin, the 

O 

density, p, is in g/cm , and the viscosity n, is in g/(cm-sec). 

The uncertainty in the lower temperature range, using the above 
formula, is estimated in reference 7 to be approximately one-half 
percent and increases to about ten percent if the same formula is 
used for temperatures up to 300 K. 
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APPEiroiX H 

Vapor Pressvire for Parahydrogen 

The vapor pressiire is the pressure, as a function of temperature, 
of a liquid in equilibrium with its own vapor. For parahydrogen, from 
reference 6, 

P 

P = 10 ° (T-29)^ + Ag(T-29)^ + A,j,(T-29)'^ (ll.l) 


vhere 


P 

o 



^ — + A 

T + A„ !)• 


T 


(H.2) 


and 


A^ = 2.000620 A^ = -50.09708 

A^ = 1.004U Aj^ = 1.748495 x lO"^ 


0 for T £ 29 K 

1.317 X 10"^ for T > 29 K 



0 for T £ 29 K 

-5.926 X 10“^ for T > 29 K 


0 for T £ 29 K 
3.913 X 10“ for T > 29 K 

The units of pressure, P, is atmospheres while the units on the 
temperature, T, is Kelvin. 
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APPENDIX I 

Specific Cooling Capacity of Parahydrogen 

The specific cooling capacity is the amount of heat which may he 
ahsorhed hy a unit mass of parahydrogen in "being warmed from storage 
conditions to the final tunnel total conditions. In this investigation 
the assumption is made that the liquid hydrogen is stored at 
atmospheric presstire. The liquid is then compressed to the operating 
pressure of the wind tunnel before being injected into the tunnel 
circuit. Assuming the compression is isentropic and assuming the 
vaporizing and warming of the injected hydrogen occurs at a constant 
pressure, the specific cooling capacity, as given in figure 52 , is 
determined by 


specific cooling capacity = I T dS 


ity = £ 


(I.l) 


= E H 
t “l 


where is the enthalpy at the tunnel total conditions and ‘ 

is the enthalpy at the liquid initial conditions. 

The specific cooling capacity for parahydrogen, in hj/kg, may 
be approximated by 

5 . 

specific cooling capacity = ^ E.T (1.2) 

i=0 ^ 


where 


T = temperature, K 



and vhere the values of 


are given hy 


B. 

1 


= E 

j=0 


B, .P'J 


(1.3) 


where 


P = pressure, atmospheres 


and the values of B, . 

10 

Bqq = 2.3362 X 10^ 

= 1.2591 X 10^ 

Bgo = -7.U035 X 10"^ 

B^q = 9.0181 X 10"^ 

= -3.3958 X 10“^ 

== U.2281 X 10“^ 


are as follows: 


‘01- 

= -2.6924 

X 10^ 

'11 

= 7.4355 X 

10“^ 

^21 

= -8.9004 

X 10“^ 

'31 

= 5.2144 X 

ir\ 

1 

0 
' H 

‘l;l 

= -1.4658 

X lO""^ 

‘51 

= 1.5789 X 

10-"" 


= -1-9693 

B^2 = T.1273 X 10"^ 

B_ = -9.68U3 X 10"^ 

C^C. 

B^g “ 6.0l7i» X 10"^ 

Bjj2 = -1.7605 X IC"® 

B^2 " 1-9350 X 10“^^ 


The values of the specific cooling capacity calculated froui 
equation 1,2 agree >rith the values calculated from equation I.l 
generally within two percent for pressures from one to ten atmospheres 
and for temperat\ires from saturation to 300 K* 



